MgN: a new promising material for spintronic applications 
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Density functional theory calculations demonstrate that rocksalt MgN is a magnetic material at 
the verge of half-metallicity, with an electronic structure robust against strong correlations and 
spin-orbit interaction. Furthermore the calculated heat of formation describes the compound as 
metastable and suggests that it can be fabricated by tuning the relative Mg and N abundance 
during growth. Intriguingly the equilibrium lattice constant is close to that of MgO, so that MgN 
is likely to form as an inclusion during the fabrication of N-doped MgO. We then speculate that 
the MgO/MgN system may represent a unique materials platform for magnetic tunnel junctions not 
incorporating any transition metals. 



The so called d° magnets are currently challenging our 
conventional understanding of magnetism. These are 
wide-gap semiconductors and oxide insulators display- 
ing magnetic properties which cannot be ascribed to the 
presence of partially filled d or f shells^. Prototypical 
examples are diamagnetic materials in their well crystal- 
lized bulk form, which present magnetic features when 
grown as defective thin films. A similar phenomenology 
is found also in oxides doped with light elements such as 
N or Unfortunately, despite the increasingly large 

number of reports of dP magnets, their typical experi- 
mental characterization is often limited to magnetometry 
with little information about the local electronic struc- 
ture. The phenomenon thus remains rather obscure. It 
is characterized by a poor degree of reproducibility and 
still lacks of a convincing theoretical framework. 

Most of the theoretical work to date is based on den- 
sity functional theory (DFT) utilizing local approxima- 
tions to the exchange and correlation functional (LDA 
or GGA). In most of the cases the LDA/GGA, probably 
correctly, associates the magnetic moment formation to 
spin-polarized holes residing on either cation vacancies 
or impurities. In addition the ground state is always pre- 
dicted ferromagnetic and rather robust. However LDA 
and GGA usually fail in describing the symmetry de- 
tails of the wave-function, in particular the polaronic 
distortion around the magnetic center. Therefore more 
advanced techniques, including strong electron correla- 
tion, appears as necessarj^jSil. When applied, these 
confirm the nature of the magnetic moment formation 
but almost always predict no long-range magnetic cou- 
pling and no ferromagnetism. MgO:N is a prototypical 
cxample^^. Quantum Monte Carlo calculations for 
the Anderson-Haldanc model indicates the possibility of 
high temperature ferromagnetism^. However DFT cal- 
culations using either the LDA+?7 or the self-interaction 
correction (SIC) scheme return impurity levels deep in 
the gap, holes trapped by polaronic distortions and no 
room temperature ferromagnetismii£i£. 

Finally another important, but yet scarcely studied, 
class of d° magnets is represented by zincblende (ZB) II- 
V or II-IV compounds such as CaP, CaAs, CaSb±2.iIiI2 
or CaC, SrC and BaC— . These are weakly covalent solids 
with either one or two holes per formula unit, and they 



are generally predicted to be half- metallic. Intriguingly, 
although at least for the II-V's the thermodynamical sta- 
ble phase has the stoichiometric Zn3P2-type structure^, 
their rocksalt (RS) phase has usually a lower energy than 
the ZB and it is predicted metastable. 

In this work we investigate the magnetic ground state 
of RS MgN, which was already predicted as half-metal 
in the ZB phased. This is an extremely interesting ma- 
terial because of its structural similarity to the parental 
MgO, widely used as tunnel barrier in magnetic tunnel 
junctions^. We will demonstrate that MgN is a half- 
metal with negligible spin-orbit interaction and negative 
formation enthalpy. This means that, although it is not 
the thermodynamical stable phase, it can be grown as a 
metastable compound in the form of either thin films or 
as nanoclusters inside MgO matrices. 

Our calculations are performed with a development 
version of the siesta codei£, which includes the atomic 
self interaction correction (ASIC)i£ and the LBA+U— 
schemes. The core electrons are treated with norm- 
conserving Troullier-Martin pseudopotentials and the va- 
lence charge density is expanded over a numerical orbital 
basis set, including double-^ and polarized functions^. 
The real space grid has an equivalent cutoff larger than 
600 Ry, and we sample a minimum of 20 x 20 x 20 fc-points 
in the Brillouin zone. The atomic coordinates are relaxed 
by conjugate gradient until the forces are smaller than 
0.01 eV/A. For the evaluation of the spin-orbit coupling 
we have performed full-potential linear augmented plane- 
waves (FLAPW) calculations^ with the scheme imple- 
mented in FLEUR code2£. We use a 12 x 12 x 12 fc-point 
mesh and a plane wave expansion with /c max = 4 a.u. 
A muffin-tin sphere of a radius 2.02 Bohr was taken for 
both Mg and N with £ max —8 as cutoff in the ^-expansion 
of the muffin orbitals. Spin-orbit is included by second 
variation method^. 

We start our analysis by briefly comparing the ground 
state of the ZB and RS structures. In table [J we report 
the equilibrium lattice constant, a, the magnetic moment 
per formula unit, [i, and the magnetic moment formation 
energy, AEm, for the two structures. AEm is defined as 
the energy difference per formula unit between the spin- 
polarized and the non-spin-polarizcd DFT solution and 
it is positive when the magnetic state has a lower en- 
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crgy. Importantly AEm is always substantially larger 
than zero, demonstrating that the magnetism is robust 
for both the structures. We also note that our results 
for the ZB phase are in excellent agreement with previ- 
ous calculations^, and that for the RS phase AEm is 
slightly smaller than that calculated for the same phase 
of SrN and CaNi^, although it is still well above room 
temperature. 
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TABLE I: Equilibrium lattice constant a (in A), magnetic 
moment per formula unit fi (in /xb) and AEm (in meV) for 
ZB and RS MgN. 

Next we discuss in details the electronic structure of RS 
MgN. The GGA bandstructure calculated at the equilib- 
rium lattice constant of 4.44 A and displayed in Fig. [TJ 
clearly reveals the half-metallic character of MgN with a 
direct gap of 5 eV at the T point in the minority spin- 
channel. This is different from the bandstructure for RS 
CaN where the gap is indirect with the conduction band 
minimum located at Xii. The differences between CaN 
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FIG. 1: GGA band structure for RS MgN calculated at the 
equilibrium lattice constant of 4.44 A: panel (a) majority spin, 
panel (c) minority. The horizontal red line denotes the po- 
sition of Fermi level. In panel (b) we report the associated 
total density of states (black solid line) and density of states 
projected over the N-2p (red dashed line) and Mg-3s (green 
shadowed area) 

and MgN can be understood by looking at the orbital pro- 
jected density of states (PDOS) presented in Fig. [ljb). 
In the case of MgN the N 3d levels are unbound and they 
only weakly mix with the Mg 3s, which form the conduc- 



tion band. In contrast the empty Ca 3d levels provide the 
dominant contribution to the conduction band of CaN, 
shifting the band minimum to the X point. Furthermore 
their hybridization with the N 2p orbitals forming the va- 
lence band is rather strong and result in the CaN valence 
band being considerably more flat than that of MgN. In 
summary the main difference between RS MgN and CaN 
is in the different role played by the cation 3c? orbitals. 
Interestingly it was suggested^ that the contribution of 
the Ca-d level in Ca pnictides is fundamental to promote 
the half metallicity. Later such a claimed was challenged 
by the statement that the magnetic moment formation in 
II-V compounds is driven by the strong atomic character 
of the anions^. MgN demonstrates that both of the ar- 
guments are partially correct. Indeed the ionicity in SrN, 
CaN and MgN is able to localize the hole on N leading 
to the formation of the magnetic moment. However the 
hybridization with the d shell of the cation increases the 
flatness of the top of valence band causing an enhance- 
ment of AEm- 

We then investigate the robustness of the half- 
metallic ground state with respect to compression, spin- 
fuctuations, choice of exchange and correlation functional 
and spin-orbit interaction. In figure[2]we report the mag- 
netic moment as a function of the lattice constant as cal- 
culated with GGA (the LDA results are rather similar). 
Importantly we note that there is a broad range of lat- 
tice constants for which MgN preserves the half-metal 
moment of 1 /is- For a smaller than 4.3 A the moment 
becomes fractional and the material turns into a standard 
ferromagnet, until it finally becomes non-spin-polarized 
for a massive compressive pressure (a < 3.0 A). Inter- 
estingly at the MgO lattice constant (4.21 A) the mag- 
netic state of MgN is still extremely close to that of a 
half-metal. This is a tantalizing feature suggesting that, 
should MgN be made in MgO, it will be an efficient spin 
polarizer/analyzer in transition metal free tunnel junc- 
tions. Finally note that for a largely expanded lattice 
parameter the moment reaches the value of 3 /ib expected 
from isolated N ions. 

Spin excitations and the possible effects arising from 
strong correlations arc analyzed next. In particular we 
calculate the total energy difference per formula unit be- 
tween the ferromagnetic and either the type-I (AE^_^) 
and type-II (AE^_^) antiferromagnetic state. These 
are both characterized by a planar ferromagnetic align- 
ment with antiparallel orientation between planes, ei- 
ther along the (001) (type-I) or (111) (type-II) direction. 
At the equilibrium lattice constant we find AE^_^ = 

52 meV (AE^_ n = 67 meV) and for AE]\_ n = 

53 meV (AE\\_ n = 68 meV) for LDA (GGA). These 
are rather similar to AEm indicating that in MgN Stoner 
and spin- wave excitations compete. 

It is then interesting to investigate how the mo- 
ment formation and the energetic of the spin-excitations 
change when electron correlation is added to the elec- 
tronic structure. This is investigated with the LDA+C7 
scheme applied to the N-2p shell as a function of U (J 
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FIG. 2: Magnetic moment n per formula unit as a function of 
the lattice constant for RS MgN as calculated with GGA. The 
inset shows a zoom at around the GGA equilibrium value of 
4.44 A. The arrows indicate the experimental MgO and the 
GGA-calculated MgN lattice constants. 



is kept at zero) . Our results are reported in table [TTJ As 
expected the addition of on-site correlation to the N-2p 
shell downshifts its center of mass effectively opening the 
band-gap in both the spin-bands. This is associated to a 
considerable strengthening of the ferromagnetic interac- 
tion, with both AEl^_+^ and AS||_^ increasing mono- 
tonically with U. In the table we also report the exchange 
split, AE CX , which is the energy difference between the 
top of the majority and minority valence bands for the 
ferromagnetic ground state at its equilibrium lattice con- 
stant. We note that also AE 0X increases with U with an 
almost perfectly linear dependence. Finally we find that 
the critical lattice constant for the half-metallicity gets 
progressively reduced as U gets larger, while the equi- 
librium lattice parameter seems to be only marginally 
affected by U. 

The promotion of the magnetic stability in RS MgN 
with increasing U can be understood by using the DFT 
version of the standard Stoner argument applied to the 
LDA+J7 functional^. In fact, in the limit of uniform 
occupancy of the strongly correlated shell (N-2p in our 
case), one expects the Stoner parameter to increase lin- 
early with U. This may, of course, be compensated by 
a reduction of the density of states at the Fermi en- 
ergy (E-p), resulting in an overall reduction of the mag- 
netic interaction. In the case of RS MgN, however, both 
the band-width and the curvature of the bands near E-p 
change little with U , since such a valence band has a pure 
N-2p character. As a consequence we find an almost per- 
fect linear relation between the exchange split AE CX and 
U , which confirms the Stoner nature of the ferromag- 
netism in MgN. This leads us to conclude that the ad- 
dition of strong correlations through an on-site Coulomb 
repulsion enhances the stability of the half-metal charac- 
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TABLE II: Critical lattice constant for the half-metallicity, 
a c , minority spin bandgap, -E gap , exchange split A£W, and 
ferromagnetic to antiferromagnetic total energy difference 
AEl-t_-t\ and Ai5f|_ Ti , as a function of the on-site Coulomb 
interaction U. 

In table HI] we also report results obtained with the 
ASIC method, which was proved to yield an accurate 
MgO bandgap and the proper polaronic distortion in 
diluted d° magnets^. The ASIC returns both a criti- 
cal lattice constant for the half-metallicity and magnetic 
energies similar to those calculated with LDA+t/ and 
U ~ 3 eV. However the bandgap and A ex are consider- 
ably larger and compatible with U = 6 eV. Although it is 
difficult to make any strong quantitative statement, be- 
cause the mixed success of ASIC in predicting exchange 
constants^, the result further confirms that the half- 
metallic state of MgN is strong with respect to strong 
electron correlation. 

Finally we investigate the effect of the spin-orbit cou- 
pling on the electronic structure of MgN. In general spin- 
orbit interaction always destroys the half-metallicity of a 
magnet by mixing the two spin channels. Our calcu- 
lations however show that the spin-orbit coupling con- 
stant between the spin and the orbital angular momen- 
tum is only of the order of 10 meV. Furthermore, the 
MgN highly symmetric cubic enviroment quenches the 
angular magnetic moment, which is only 7 x 10 -5 /ib for 
Mg and 9.8x10 _4 /zb for N. It is then not surprising that 
the spin polarization remains 100 % even after the spin- 
orbit interaction is accounted for, i.e. the half-metallicity 
is preserved. 

After having discussed the electronic properties of 
MgN we now turn our attention to analyzing the struc- 
tural stability of the RS structure. First we compare the 
total energy per cell of volume of four different crystal 
structures, likely candidates for MgN, namely wurzitc, 
ZB, NiAs-typc (space group PG^/mmc) and RS. As 
shown in figure[3]thc RS structure presents both the lower 
total energy and smaller equilibrium volume of all the 
crystals investigated, indicating that it is the most ther- 
modynamic favorable. This however still does not mean 
that MgN can be made. In fact one has to compare the 
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FIG. 3: (Color on line) Total energy as a function of the vol- 
ume per formula unit for MgN with wurtzite (WZ), zincblende 
(ZB), NiAs and rock-salt (RS) crystal structure. 



stability of RS MgN with that of its stechiometric stable 
form Mg3N2, i.e. one has to demonstrate that the heat 
of formation for RS MgN, Ai/MgN, is comparable to that 
of Mg 3 N 2 , Ai? Mg3N2 . 

Mg 3 N 2 is a non-magnetic semiconductor with an ex- 
perimental band gap of 2.8 eV— and with the tetragonal 
bixbyite crystal structure. The experimental lattice con- 
stant is 9.953 AM while our calculations return values of 
9.85 A (LDA) and 10.01 A (GGA). These are both within 
1% of the experimental reported values. The heats of for- 
mation can be written as 

Ai?Mg 3 N 2 = 1/5 [1/8 i?Mg 3 N 2 — (3/2 i?Mg + £n 2 )] i (1) 

Ai/MgN = 1/2 [i?MgN — 1/2 [Eu s + -E-n 2 )] j (2) 



where E refers to the total energy and Mg is assumed 
in its hep metallic phase. We calculate A77Mg 3 N 2 = 
— 1.44 eV/atom (Ai/Mg 3 N 2 = —1.16 cV/atom) and 
AiJ M gN = -0.91 eV/atom (Ai? M gN = -0.59 eV/atom) 
with LDA (GGA). Although, as expected Mg 3 N 2 has a 
lower heat of formation and forms at equilibrium, RS 
MgN also has a negative AH. This means that it is 
a thermodynamic stable compound, and hence that RS 
MgN is metastable. We then speculate that the relative 
abundance of N and Mg during the growth process, in 
addition to the choice of substrate, can drive the forma- 
tion of either Mg 3 N 2 or MgN. Intriguingly MgO appears 
as an ideal substrate for the growth of MgN but metals 
with bec lattice structure are equally tantalizing. Among 
them we mention V, Mo and W. 



In conclusion we investigated the electronic and struc- 
tural properties of rocksalt MgN. This was found to be 
at the verge of half-metallicity in LDA and GGA and 
completely half-metallic as soon as some electron corre- 
lation is included. The half-metallicity is then robust 
with compression and spin-orbit interaction. Further- 
more we have then demonstrated that the RS phase is 
the most stable among the various possible MgN phases. 
Its heat of formation is negative, but smaller than that 
of Mg 3 N 2 , suggesting that MgN can be possibly made by 
tuning the relative Mg and N abundances during growth 
and the substrate crystal structure. 
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